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NOL  HYPERSONIC  TUNNEL  NO.  4 RESULTS  VII: 
EXPERIMENTAL  INVESTIGATION  OF  TURBULENT 
BOUNDARY  LAYERS  IN  HYPERSONIC  FLOW 


Prepared  by: 

R.  Kenneth  Lobb,  Eva  M.  Winkler,  and  Jerome  Persh 


ABSTRACT:  Naturally  turbulent  boundary  layers  on  the  wall 
of  a wedge-typo  water-cooled  nozzle  in  the  NOL  12  x 12  cm 
Hypersonic  Tunnel  No.  4 have  been  investigated  at  Mach 
numbers  of  5.0  to  8.2  with  and  without  steady  state  heat 
transfer  to  the  surface.  The  Reynolds  number  based  on 
boundary  layer  momentum  thickness  v/as  varied  from  5,000  to 
13,000.  Measurements  of  Pitot  and  static  pressures,  total 
and  wall  temperatures,  and  rates  of  heat  transfer  made  it 
possible  to  compute  velocity  profiles,  temperature  profiles, 
and  boundary  layer  parameters  without  resorting  to  any 
assumptions.  The  turbulent  portion  of  the  boundary  layer 
velocity  profile  was  found  to  differ  from  the  incompressille 
flow  logarithmic  law  by  an  amount  that  depends  on  the  heat 
transfer  and  Mach  number.  The  data  for  the  outer  turbulent 
portion  for  any  one  Mach  number  fall  on  a single  curve  if 
plotted  in  a particular  nondimensionai  coordinate  system 
which  is  based  on  local  properties  in  the  boundary  layer. 

The  velocity  profile  in  the  laminar  sublayer  is  linear.  The 
thermal  sublayer  was  found  in  all  cases  to  be  larger  than 
the  velocity  sublayer.  Local  skin  friction  coefficients 
for  zero  heat  transfer  as  calculated  from  the  velocity  gradi' 
ents  at  the  wall  arc  consistent  with  the  results  of  other 
experimenters  at  lower  Mach  numbers  for  the  same  Reynolds 
number.  Furthermore  it  is  demonstrated  that  the  Reynolds 
analogy  between  skin  friction  and  heat  transfer  is  valid  at 
aiach  numbers  up  to  8.2. 
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This  is  the  seventh  NAVORD  Report  on  an  investigation 
carried  out  in  the  continuous  NOL  12  x 12  cm  Hypersonic 
Tunnel  No.  4.  The  titles  of  the  previous  NAVORD’s  describing 
results  from  the  tunnel  are : 

I.  Air  Liquefaction.  NAVORD  1742,  4 January  1951 

II.  Diffuser  Investigation.  NAVORD  2376,  5 May  1952 

I I I.  Diffuser  Investigation  with  Models  and  Support. 
NAVORD  2435,  1 July  1952 

IV.  High  Supply  Temperature  Measurement  and  Control. 
NAVORD  2574,  8 October  1952 

V.  Experimental  and  Theoretical  Investigation  of  a 
Cooled  Hypersonic  Wedge  Nozzle.  NAVORD  2701, 

13  April  1953 

VI.  Experimental  and  Theoretical  Investigation  of  the 
Boundary  Layer  and  Heat  Transfer  Characteristics 
of  a Cooled  Hypersonic  Wedge  Nozzle  at  a Mach 
Number  of  5.5.  NAVORD  3757,  8 July  1954 

The  present  NAVORD  (Results  VII)  presents  experimental 
investigations  and  discussions  of  turbulent  boundary  layers 
in  hypersonic  flow. 

Knowledge  of  the  related  effect  of  skin  friction  and  aero- 
dynamic heating  at  high  Mach  numbers  is  needed  by  the 
designers  of  hypersonic  vehicles  because,  for  slender 
missiles,  skin  friction  represents  the  major  part  of  the 
drag  and  extreme  surface  temperatures  produced  by  frictional 
heating  may  very  well  be  a limiting  design  factor. 

A portion  of  the  results  contained  in  this  NAVORD  were 
presented  at  the  22nd  Annual  Meeting  of  the  Institute  of 
Aeronautical  Sciences  and  also  at  the  Bureau  of  Ordnance 
Committee  on  Acroballistics  Symposium  in  October  1954.  The 
present  report  contains  additional  results  and  more  detailed 
analysis  of  the  data  a.s  well  as  a complete  tabulation  of  the 
experimental  results. 

This  work  was  jointly  sponsored  by  the  U.  S.  Naval  Bureau  of 
Ordnance  an*’  the;  U.  S.  Air  Force.  It  was  carried  out  under 
Tasks  N0L-M9a-108-l-54,  NOL-M9a-1.33-l-55 . and 
NOL-M9a-133-5-55. 

The  authors  are  indebted  to  Dr.  R.  E.  Wilson  for  many  stimu- 
lating discussions  during  the  course  of  the  investigations. 


ii 


NAVORD  Report  3880 


The  numerical  evaluation  of  the  boundary  layer  surveys  was 
done  on  a card-programmed  calculator  of  the  Applied  Mathe- 
matics Division.  The  preparation  of  the  data  for  the  card- 
programmed  calculator  was  greatly  assisted  by  Dr.  fi.  K.  Blum. 
The  cooperation  of  Messrs.  L.  L.  Liccini  and  R.  Garren,  Jr., 
who  participated  during  the  tests,  is  acknowledged. 


JOHN  T.  H.AYWARD 
Captain,  USN 
Commander 


H.  H.  KURZWEG,  Chief 
Aeroballistic  Research  Department 
By  direction 
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SYUBOLS 


V/S.rC*^  ;d«vJLli  JLX^CCJLOi*  A V w«I&wL<  A A ^ w **o  La  Cm<u 

2T 

conditions  


w 


Pod  '*cd 


Cji  incompressible  local  skin  friction  coefficient  for 
zero  heat  transfer  bt'iseu  on  free-streara  conditions 

Cp  specific  heat  at  constant  pressure 

H boundary  layer  shape  parameter  /^*/9 

k ti'.ermal  conductivity 

M Mac.',  -imber 

Nu  Nusselt  . ' or  for  stagnation  temperature  probe 

Pr  Prandtl  numbet 

stagnation  pressure 

p^  Pitot  pressure 


p static  pressure 

II  Gas  constant 

Re  Reynolds  numl'or  based  on  iree-stream  conditions 

r recovery  factor 

St  Stanton  number  (equation  (3)) 

local  stagnation  temperature 

stagnation  temperature  as  measured  by  a stagnation 
temperature  probe 

T local  static  temperature 

u velocity 

u'  velocity  component  perpendicular  to  wall 

I 

friction  velocity  {.V^/p)^ 
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u"*"  velocity  parameter  u/u^  ( p based  on  wall  properties) 

u*  velocity  parameter  u/u^  ( p based  on  local  properties) 

y distance  perpendicular  to  wall 

y"^  wall  distance  parameter  y /**'  (p  and  based  on 
wall  properties) 

y*  wall  distance  parameter  y Uy  ( p and  X based  on 
local  properties) 

^ ratio  of  specific  heats 

(f  total  boundary  layer  thickness 

Q * boundary  layer  displacement  thickness 

<r 

0 momentur.i  tliicknoss  / - — fl  - rr — ^]dy 

^ JPOO^OD  t 


P 

r 


viscosity 

kinematic  viscosity 

density 

shear  stress 


Subscripts 

e equilibrium  wall  temperature  for  zero  heat  transfer 
n physical  properties  of  gas 

L edge  of  laminar  sublayer 

p temperature  probe 

L physical  properties  of  thermocouple  wire 

T temperature  profile 

u velocity  profile 

w values  based  on  wall  conditions 

0 values  based  on  momentum  thickness 
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values  bar 'd  on  free-stream  conditions  outside 
boundary  layer 

conditions  behind  shock  in  front  of  probe 
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■-(OL  HYPERSONIC  TUNNEL  NO.  4 RESULTS  VII; 
EXP’iRIMENTAL  INVESTIGATION  OF  TURBULENT 
BOUNDARY  LAYERS  IN  HYPERSONIC  FLOW- 


INTRODUCTION 

1 . Boundary  layer  investigations  at  hypersonic  speeds  are 
of  immediate  practical  interest  since  friction  drag  and 
heat  transfer  data  are  needed  by  the  designers  of  hypersonic 
missiles . 

2.  From  the  theoretical  viewpoint,  various  treatments  of  the 
problem  are  available  for  the  case  of  laminar  boundary  layer 
flows  as  well  as  for  turbulent  boundary  layers.  The  experi- 
mental data  are,  however,  limited,  and  in  general  cover  the 
Mach  number  range  up  to  about  5.  The  lack  of  experimental 
confirmation  for  laminar  boundary  layer  theory  is  not 
regarded  as  serious  since  the  theory  is  considered  to  be  more 
or  less  exact.  Most  of  the  analyses  of  compressible  turbu- 
lent boundary  layer  flows,  however,  are  based  on  experimental 
results  obtained  in  wholly  incompressible  flow.  Tne  uncer- 
tainties in  those  analyses  have  led  to  large  discrepancies 

in  the  prediction  of  skin  friction  and  heat  transfer.  Even 
those  theories  which  agree  with  each  other  and  with  the  exist- 
ing experimental  data  within  5 percent  up  to  Mach  numbers  of 
5 differ  greatly  at  hypersonic  Mach  numbers  (reference  (a)). 

3.  The  purpose  of  the  present  investigation  was  to  extono 
the  Mach  number  range  of  available  data  and  also  to  attempt 
to  provide  a deeper  insight  into  the  characteristics  of  a 
turbulent  boundary  layer  in  compressible  flows. 


Experimental  Equipment  and  Techniques 

4.  The  boundary  layer  surveys  have  been  conducted  in  the 
NOl  12  X 12  cm  Hypersonic  Tunnel  No.  4,  which  is  described  in 
references  (b) , (c) , and  (d).  This  tunnel  operates  contin- 
uously in  the  Mach  number  range  from  5 to  10  free  of  air 
condensation  effects.  Supply  temperatures  from  300  to  800°  K 
and  supply  pressures  from  1 to  15  atmospheres  are  available. 

A two-d imr ns i onal  watcr-coolcu  wedge  nozzle  expands  the  air 
and  can  be  adjusted  for  any  Mach  number  simply  by  changing 
the  tliroat  oponini,  area. 

5.  The  investigations  were  made  for  a ran  ,c  of  Mach  numbers 
varying  frem  5 to  about  S.  For  each  Mach  number  several 
surveys  were  mauc  diffcrin..  ^n  the  Reynolds  number  and  the 
rate  of  heat  transfer  to  the  wall.  The  surveys  were  made 

at  the  conti.r  of  one  nozzle  wall,  appro.ximatoly  four  inciies 
upslro.'im  (jf  lh<'  nozzle  exit  plane.  Fi-uro  1.  Measurements 
at  this  pus-l.o:i  arc  unaffected  Ijy  tlic  junction  of  tiie  nozzle 
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e,iu'  and  the  first  diffuser  plate  (reference  (e)).  The  rela- 
tively thick  (approximately  25  mia)  turbulent  boundary  layer 
on  the  nozzle  wall  minimizes  probe  positioning  errors  and 
thus  facilitates  accurate  determination  of  the  profile  shapes. 
Such  a boundary  layer,  however,  is  sub.iected  to  a slight 
frcc-stream  pressure  gradient  in  the  neighborhood  of  the 
survey  s La  lion  because  the  flow  is  radial.  (The  corresponding 
Mach  numi'cr  rise  is  about  3 percent  per  tunnel  caliber  at  the 
survey  plane  for  a Mach  number  of  5 and  decreases  with 
increasing  Mach  numher.)  it  is  felt  that  the  effect  of  the 
pressure  gradient  on  the  profiles  is  small.  The  wall  tempera- 
ture, I the  other  hand,  is  maintained  practically  constant 
near  room  temporaturo  over  the  entire  length  of  the  nozzle, 
except  at  the  throat  where  the  surface  temperature  approaches 
the  recovery  temperature.  In  all  tests  the  transition  from 
laminar  to  turbulent  boundary  layer  occurred  slightly  down- 
stream of  the  nozzle  throat*.  The  exact  effect  of  the  history 
on  the  local  boundary  layer  charact-eristics  is,  of  course, 
not  fully  known  and  the  results  must  be  interpreted  with  this 
fact  in  mind. 


G.  Reynolds  numbers  and  the  rate  of  heat  transfer  to  the 
wall  are  controlled  by  the  supply  pressure  and  the  supply 
temperature.  In  all  tests  the  Reynolds  numbers  based  on 
boundary  layer  thickness  .are  roughly  two  orders  of  magni- 
tude greater  tiian  needed  to  make  slip  flow  effects  at  the 
surface  negligible  (reference  (f)).  Since  the  wall  tempera- 
ture is  always  maintained  near  room  temperature,  the  lowest 
rate  of  heat  transfer  to  the  water-cooled  nozzle  that  could 
;.(■  realized  corresponds  to  a supply  temperature  which  is  just 
hi;,!!  oMuugii  to  avoid  air  condensation  in  the  test  section. 

7.  For  oacl;  survey,  Pitot  and  static  pressures,  stagnation 
.nnd  wall  temperatures,  and  the  wail  temperature  gradient 
poriiondicular  to  the  nozzle  surface  are  recorded.  Static 
pressures  are  measured  by  a 0.G4  mm  diameter  orifice  in  tlie 
wall  and  by  a static  probe  in  the  free  stream  just  outside 
the  £*dge  of  the  boundary  layer.  The  static  probe  is  an  8° 
cone  cylinder  with  orifices  located  17  diameters  aft  of  the 
slioulcier,  Fi,gure  2.  It  is  mounted  from  the  side  wall  with 
its  axis  parallel  to  the  flow.  Since  agroemont  between 
wall-  .Tnd  free-stroam  static  pressures  is  within  1 percent, 
a constant  static  pressure  has  been  assumed  to  exist  through 
the  boundary  layer.  The  measurements  arc  made  with  silicone 
oil  manometers  of±2  microns  measuring  accuracy. 


‘Surface  probe  tests  made  subsequent  to  the  results  reported 
}!'.  reference  (e)  indicated  that  the  boundary  layer  profile 
reported  in  this  reference  was  close  to  the  transition  region. 
It  was  definitely  established,  from  the  results  of  the  sur- 
face probe  tests,  that  all  boundary  layei  profiles  presented 
herein  were  measured  in  the  fully  turbulent  region. 
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8.  The  Pitot  pressure  is  surveyed  from  wall  to  free  stream 
with  a flattened  hypodermic  tube  of  0.125  mm  half  height, 
Figure  2.  The  opening  is  large  enough  to  avoid  errors  due 
to  slip  flow  effects  on  the  pressure  measurements  in  the 
region  of  low  Reynolds  number  in  the  boundary  layer  near  the 
wall  (reference  (g)).  The  probe  is  mounted  in  a micro- 
traverse  mechanism  which  also  accommodates  the  connection 

to  the  pressure  gauge.  The  position  of  the  centerline  of 
the  probe  is  measured  relative  to  the  position  of  electrical 
contact  between  probe  and  wall  with  an  accuracy  of  ±0.025  mm. 
Impact  pressures  above  20  mm  Hg  are  measured  with  a precision 
mercury  manometer  of  10.1  mm  measuring  accuracy.  For  the 
lower  pressiires  a silicone  oil  manometer  is  used. 

9.  A stagnation  temperature  probe  with  a single  platinum- 
coated  silica  shield  is  used  for  the  temperature  surveys 
of  the  boundary  layer,  Figure  2 (reference  (h)).  For  the 
neasurements  close  to  the  wall,  a flattened  probe  with  a 
half  height  of  0.48  mm  is  used.  The  temperature  recovery 
factor  of  tho  probes  reaches  a value  of  0.998  for  large 
Reynolds  numbers.  The  variation  of  the  temperature  recovery 
factor  with  flow  parameters  for  each  probe  is  described  by 

a single  calibration  curve  by  relating  the  calibration  data 
to  the  flow  conditions  inside  the  probe.  Reference  to  this 
curve  makes  it  possible  to  determine  eccurately  total  and 
static  temperatures  throughout  the  boundary  layer  from 
measured  temperatures  and  pressures.  The  e.m.f.  output  of 
the  temperature  probe  is  recorded  on  Brown  temperature 
recorders.  Tho  temperature  measurements  are  accurate  to 
t0.2°  0. 


10.  LiOcal  values  of  the  heat  transfer  to  tho  nozzle  wall  and 
nozzle  surface  temperatures  at  the  boundary  layer  survey 
station  are  obtained  from  temperature  measurements  in  tho 
nozzle  wall.  These  measurements  are  made  with  four  thermo- 
couples imbedded  in  the  nozzle  wall  at  various  distances 
from  the  surface  (reference  (d)).  Depending  upon  the 
operational  conditions  of  the  tunnel  and  the  rate  of  coolant 
flow  to  the  nozzle,  equilibrium  readings  are  reached  after 
10  to  20  minutes  of  tunnel  operation.  They  show  that  the 
temperature  drops  linearly  from  the  nozzle  surface.  Previous 
investigations  indicate  that  longitudinal  and  lateral  tempera- 
ture gradients  at  the  nozzle  surface  can  be  neglected 
(reference  (d)).  The  wall  temperatures  are  measured  to  within 
10.01°  C with  a Leeds  and  Northrup  K-2  type  potentiometer. 


Data  Reduction 

11.  Tho  experimental  data  needed  to  evaluate  tho  boundary 
layer  profiles  are: 

a.  Pitot  pressure,  p^,  vs.  distance  from  wall 
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b.  Static  pressure,  p,  taken  constant 

c.  Ifeasured  total  temperature,  Tj^,  vs.  distance  from 
wall 

d.  Variation  of  probe  temperature  recovery  factor  with 
flow  parameters 

Since  p^  and  are  not  always  measiured  at  identical 

positions  from  the  wall,  both  sets  of  data  are  taken  at 
sufficiently  close  intervals  so  that  the  temperatures  corre- 
sponding to  impact  pressure  readings  at  a particular  position 
can  be  interpolated  from  a curve  faired  through  the  tempera- 
ture data. 

12 . The  numerical  evaluation  of  the  data  is  done  on  the  HOL 
Card-Programmed  Calculator.  The  Rayleigh  formula  is  used  to 
compute  the  Mach  number  profile  from  p^  and  p (reference 

(i)).  To  evaluate  the  stagnation  temperatiure , T^,  from  the 

measured  T^  values  requires  reference  to  the  calibration 

curve  of  the  probe  (reference  (h)).  This  curve  gives  the 
recovery  factor  of  the  probe  in  terms  of  the  measured  total 
temperature  and  the  pressure,  P2,  behind  the  bow  shock  in 

front  of  the  probe 


7/4 

f(P2/RT('^’ 


CD 


13.  To  calculate  the  friction  velocity,  u<p  , for  the  u*, 
y*  representation  (reference  (J))  of  the  boundary  layer 
velocity  profiles,  the  wail  shear  stress,  V is  calculated 

either  from  the  slope  of  the  velocity  profile  in  the  laminar 
sublayer 


(2) 


or  from  the  temperature  measurements  in  the  nozzle  wall , 
using  the  following  equations: 


St 


uo 


(AT/Ay)k  

” "^w^  °p  “cD 


(3) 


4 


NAVORD  Rei>->rt  3880 


and  the  Reynolds  analogy  (reference  (k)) 


-r 


St 


00 


tv.2/3  . 2 

^ Poa  “cc 


(4) 


Values  of  21T/ Ay  are  obtained  from  the  measured  wall 
temperature  gradient  and  values  of  Tp  were  calculated  tising 
1/3 

r = Pr  with  Rr  • 0.72. 


RESULTS 


Boundary  Layer  Profiles 

14.  Naturally  turbulent  boundary  layers  were  surveyed  at 
free-streaa  Mach  numbers  of  5.0  to  8.2.  These  surveys  were 
made  at  different  Reynolds  numbers  Reg  and  various  rates 

of  heat  transfer  to  the  nozzle  wall,  as  listed  in  the 
following  table: 


“00 

Re^ 

(Te  - T^)/Te 

4.93 

5,350 

0 

5.01 

6,480 

0.223 

5.03 

7,950 

0.373 

5.06 

7,370 

0.420 

5.75 

11,600 

0.108 

5.79 

12,400 

0.239 

5.82 

11,400 

0.379 

6.83 

8,550 

0.325 

6.83 

12,640 

0.443 

6.78 

8,400 

0.437 

6.78 

7,900 

0.500 

7.67 

8,130 

0.487 

8.18 

9,540 

0.496 

A tabulation  of  the  complete  data  is  given  in  Table  I 
(page  22)  and  Appendix  A.  Typical  plots  of  the  basic  data 
are  shown  in  Figiu'es  3 through  6.  The  data  points  close  to 
the  wall  that  v/ere  influenced  by  the  presence  of  the  wall 
were  omitted  from  all  figures  and  tables. 
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15.  For  most  of  the  curves,  a dimensional  distance  has  been 
selected  for  the  aoscissa  in  order  to  show  more  clearly  the 
physical  differences  between  the  profiles  at  different  heat 
transfer  rates.  For  the  same  reason  some  of  the  curves  and 
data  points  have  ueen  omitted  from  Figures  3 through  8,  Shape 
differences  are  pronounced  only  in  the  temperatiire  profiles 
of  the  boundary  layer,  Figures  4 and  5.  The  static  tempera- 
ture profiles,  Figure  5,  obtained  with  heat  transfer  have  a 
temperature  maximum  close  to  the  wall.  The  temperature  at 
the  maximum  point  or  those  still  closer  to  the  wall  could 
not  be  evaluated  from  measured  air  temperatures  because  of 
the  physical  size  of  the  tempera t\ire  probes.  In  the  cases 
where  heat  transfer  data  have  been  measvired  in  the  nozzle 
wall  (see  Table  I),  the  slope  of  the  temperature  curve 
immediately  at  the  wall  has  been  deduced  from  the  tempera- 
ture gradient  in  the  nozzle  wall 


(5) 


(Values  of  the  thermal  conductivity  of  air  have  been  taken 
from  reference  (a)).  The  curve  through  the  measured  data 
and  the  computed  slope  have  been  joined  to  give  a smooth 
shape  to  the  static  temperature  curves.  (Since  only  the 
square  root  of  the  temperatxure  enters  into  the  computation 
of  the  velocity,  the  velocity  profile  close  to  the  wall  was 
found  to  be  insensitive  to  errors  in  T due  to  incorrect 
interpolation  between  the  last  measured  air  temperature  and 
the  wall  temperature.) 

16.  The  velocity  profiles  at  any  one  Mach  munber,  Figure  6, 
are  similar  in  the  turbulent  outer  part.  A sharp  change  in 
slope,  near  the  wall,  roughly  specifies  the  extent  of  the 
laminar  sublayer.  Within  this  sublayer  the  velocity  varies 
linearly  with  distance;  the  slope,  however,  increases  with 
increasing  heat  transfer. 

17.  It  should  be  noted  here  that  velocities  deduced  from 
Pitot  measurements  in  the  lamimar  sublayer  may  bo  slightly 
high  due  to  the  effect  of  velocity  fluctuations  on  the  Pitot 
pressure  readings  close  to  the  wall.  Since  the  fluctuating 
velocity  component,  u is  always  positive,  the  Pitot  tube 
measures  total  pressures  that  are  always  somewhat  greater 

than  the  effective  total  pressurv;  by  the  amount  J p li . 

In  Figure  7 the  variation  of  (u’  /u  with  nondimensional 
distance,  y'*‘,  from  wall  is  shown  for  the  incompressible  data 
of  references  (ra) , (n),  and  (o).  Although  there  is  consid- 
erable scatter,  the  curve  was  drawn  by  giving  the  greatest 
consideration  to  the  trend  indic.nted  by  the  larger  bulk  of 
data.  On  the  assumption  that  the  incompressible  results  of 
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Figure  7 are  applicable  to  the  present  data,  the  laminar 
sublayer  velocity  profile  data  can  be  adjusted  by  the  method 
given  in  reference  <m) . This  adjustment  indicates  that  the 
measured  velocities  at  the  edge  of  the  laminar  sublayer  may 
be  higher  than  the  true  velocity  by  a maximum  amount  of 
3 percent 


Nondimcnsional  Velocity  Profile  Representation 

18.  Turbulent  boundary  layer  velocity  profiles  arc  very 
often  presented  in  nondiraensional  coordinates  which  may  be 
computed  either  on  the  basis  of  wall  properties,  u and 
y'**,  or  local  properties  at  each  point  in  the  boundary  layer, 
u*  and  y*  (sec  list  of  symbols) . Figures  8 to  11  show 

all  the  measured  velocity  profiles  in  the  u"**,  y*  coordinate 
system.  A pronounced  upward  shift  of  the  turbulent  portion 
of  the  cvirves  with  increasing  heat  transfer  is  apparent  for 
each  Mach  number.  A similar  shift  has  been  observed  by 
Deissler  (reference  (p)}  in  his  experiments  with  subsonic 
compressible  turbulent  boundary  layers  in  pipes.  The  position 
of  the  three  upper  curves  of  Figure  10  relative  to  each  other 
indicates  that  a Reynolds  number  effect  is  probably  super- 
imposed. According  to  the  respective  values  of  the  heat 
transfer  parameters  for  these  surveys,  a different  spacing 
in  the  txirbulent  portion  would  be  e):pected. 

19.  In  addition  to  this  Reynolds  number  effect  the  data  in 
the  outer  turbulent  portion  of  the  profiles  (figures  8 to  il) 
apparently  disperse  with  Mach  number  and  heat  transfer 
parameter,  or  some  combination  thereof.  In  an  effort  to 
determine  a parameter  which  would  eliminate  the  dispersing 
effect  of  at  least  one  of  the  quantities  mentioned  above, 
the  values  of  u»  and  y*  wore  evaluated  for  each  point 

in  the  boundary  layer  for  each  of  the  measured  profiles.  A 
plot  of  these  results  indicated  a tendency  toward  a single 
curve  for  each  Mach  number,  with  the  dispersing  effect  of 
the  different  heat  transfer  rates  removed,  but  not  that  of 
the  Reynolds  number.  It  was  found  that  this  dispersion  of 
the  data  could  be  removed  by  dividing  the  y*  coordinate 
by  the  respective  values  of  the  shape  parameter  (H) . The 
data  for  any  one  Mach  number  then  converge  to  a single 
curve,  (Figure  12).  With  increasing  Mach  number  the  curves 
fan  out  in  an  upward  direction.  That  this  displacement  is 
systematic  is  demonstrated  in  Figure  13  which  shows  a cross 
plot  of  the  data  of  Figure  12.  The  naturally  turbulent 
boundary  layer  data  of  reference  (t),  as  well  as  the  M " 0 
results  are  included  in  Figure  12. 

20.  The  data  of  Figures  8 to  11  may  be  used  to  examine  some 
of  the  fundamental  assumptions  of  several  of  the  theoretical 
treatments  of  compressible  turbulent  boundary  layers.  It  is 
implied  in  these  analyses  (references  (q) , (r),  (s)  that  the 
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edge  of  tho  laminar  sublayer  occurs  at  a fixed  value  of 
u'*’  - y^,  and  that  the  slope  of  tho  turbulent  portion  of  tho 
u+i  y*  curve  is  constant,  regardless  of  Mach  number  or  boat 
transfer  conditions.  It  should  be  noted  that  the  results 
of  Figures  8 to  11  show  that  these  assujnptions  are  not 
generally  valid. 


Velocity  and  Thermal  Sublayers 

21.  It  is  apparent  from  the  curves  of  Figures  £ through  11 
that  the  u'*',  y*  velocity  profile  representation  only  roughly 
.specifies  the  extent  of  the  laminar  sublayer.  The  edge  of 
t.he  velocity  sublayer  can  be  accurately  determined  from  a 
plot  of  the  coordinates  u/uoo  against  y/<^  on  logarithmic 
paper.  On  such  a plot  it  is  found  that  the  outer  turbulent 
position  of  the  boundary  layer  may  be  well  fitted  with  a 
straight  line,  the  slope  of  which  is  the  exponent  in  the 
power  profile  representation  of  the  turbulent  boundary  layer 
profiles 


(6) 


In  the  laminar  sublayer,  a straight  line  of  unit  slope  may 
be  drawn  quite  accurately  through  the  data.  The  point  of 
intersection  of  those  two  straight  lines  is  defined  as  the 
edge  uf  the  laminar  sublayer.  Figi’re  1-!  demonstrates  the 
results  of  this  procedure  applied  to  the  data  of  several 
typical  profiles.  The  intersection  points  Indicated  are 
fairly  insensitive  to  the  slope  of  the  line  faired  through 
the  experimental  data  in  the  turbulent  portion  of  Lho  veloc- 
ity profile, and  therefore  the  scatter  of  tho  experimental 
drta  does  not  seriously  affect  the  determination  of  this 
point.  A similar  procedure  may  be  used  to  detormiae  the 
edge  of  the  thermal  sublayer.  The  thicknesses  of  the  veloc- 
ity and  thermal  sublayers  have  been  determined  for  all  of 
tho  measured  profiles  and  the  values  of  u,/u_  , (cT./eT)  , 

(cT  /</)  , and  n are  tabulated  in  Tabl®  I (page  22). 

1j  t 

22.  A comparison  between  the  values  of  {<f',/d')  and 

u 

Table  I,  shows  that  the  thermal  sublayer  is  in 

all  cases  thicker  than  the  velocity  sublayer.  According  to 
the  analysis  of  reference  (v)  the  relative  thicknesses  of 
the  thermal  and  velocity  sublayers  are  related  to  the  turbu- 
lent Prandtl  num/cr.  For  tho  present  case  of  thermal  sub- 
layers thicker  than  tho  velocity  sublayer  Reichardt’s  tnoory 
predicts  a turbulent  Prandtl  number  larger  than  the  molecular 
Ihandtl  number.  Preliminary  evaluation  of  this  quantity  fov 
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the  zero  heat  transfer  case  of  M = 4.93  and  a molecular 
Prnndtl  number  of  0.72  using  the  equations  given  by  Rubesiu 
(reference  (w))  yielded  a value  for  the  turbulent  Irandtl 
number  of  the  order  of  0.9. 


Discussion  of  Skin  Friction  Data 

23 . Because  many  of  the  theoretical  treatments  of  compress- 
ible turbulent  boundary  layers  involve  properties  of  the 
laminar  sublayer,  emphasis  was  placed  on  determining  the 
velocity  and  temperature  distributions  in  this  region. 
Furthermore  these  data  facilitate  the  accurate  determination 
of  surface  shear  stress  from  the  velocity  profile  slope  near 
the  wall . 

24.  For  most  of  the  measured  profiles,  the  surface  shear 
stress  was  calculated  from  the  slope  of  the  velocity  profile 
in  the  laminar  sublayer  and  the  heat  transfer  data  using 
Reynolds  analogy,  as  described  previously.  Values  of  c^ 
and  Stgjj  obtained  for  all  the  profiles  arc  tabulated  in 
Table  I.  In  general  the  values  of  c^  determined  from  the 

velocity  profile  data  are  up  to  5 percent  larger  than  those 
determined  from  the  heat  transfer  measurements.  A possible 
explanation  for  this  discrepancy  uay  be  given  by  considering 
the  previously  mentioned  effect  of  the  velocity  fluctuations 
on  the  Pitot  pressure  data.  Assuming  that  the  incompressible 
results  of  Figure  7 can  bo  applied,  the  laminar  sublayer 
velocity  profile  data  were  corrected  by  the  method  indicated 
in  reference  (m).  Values  of  wall  shear  stress  were  then 
computed  from  the  adjusted  velocity  profiles.  This  procedure 
brings  the  shear  stress  values  from  the  heat  transfer  data 
and  velocity  profile  data  in  closer  agreement  (.see  Table  I). 
The  two  values  are  well  within  the  experimental  accuracy  of 
the  data. 

25.  The  good  agreement  between  skin  friction  coefficients 
obtained  from  the  velocity  profile  data  and  tl>ose  determined 
from  the  measured  heat  transfer,  using  Reynolds  analogy 
demonstrates  the  applicability  of  Reynolds  analogy  at  hyper- 
sonic Mach  numbers.  That  the  Reynolds  analogy  is  applicable 
for  supersonic  flow  up  to  a i'ach  number  of  about  3.2  has  been 
shown  by  Seiff  (reference  (x)).  The  present  data,  plotted 
together  with  those  of  reference  (x),  show  that  at  hypersonic 
Mach  numbers  the  Reynolds  analogy  maintains  the  same  form 
found  to  be  accurate  for  lower  supersonic  speeds.  The  present 
data  arc  especially  convincin,:  evidence  because  the  Stanton 
number  and  the  skin  friction  <'.oof f icient  wore  ''etermincd 

s imul ianoously  by  twu  independent  experimental  methods  whereas 
the  results  of  reference  (.\)  were  calculated  by  "interpolating 
the  skin  frielion  ibti.i  Id  i-ijual  cemi  i i i (in-;  (if  Macli  number, 
wall  to  free  stream  temperature  ratio,  and  Reynolds  number. ■ 
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26.  In  the  fiixures  which  follow,  the  experimental  values 
of  Cj.  used  are  those  determined  by  the  velocity-profile 

slope  technique.  These  values  of  Cf  were  used,  not  in 

preference  to  those  obtained  from  the  heat  transfer  measure- 
ments , hut  only  because  there  arc  some  cases  for  which  heat 
transfer  measurements  are  not  available  and  it  was  desired 
to  preserve  consistency  in  the  plotted  results. 

27.  The  values  of  c^^  used  to  form  the  ratio  c^/c^^ 

were  calculated  from  the  Karsan-Schoenherr  equation  for 
the  same  Re^  value  for  zero  heat  transfer 


c 


f i 


0.586 

(logjQ2  ReQ)(logjQ2  Re^O.868) 


(7) 


The  Reynolds  number  based  on  momentum  thickness  was  used 
instead  of  the  customary  Reynolds  number  based  on  the 
distance  from  the  leading  edge.  A .’''—'olds  number  based 
on  a boundary  layer  parameter  was  selected  because  of  the 
arbitrariness  inherent  in  calculating  an  effective  leading 
edge  for  boundary  layer  measurements  on  a wind  tunnel  wall. 

28.  The  variation  of  with  heat  transfer  parumeter 

(Tg  - T^)Tg  is  shown  in  Figure  15  for  all  of  the  data  of 

the  present  investigation.  Also  shown  on  this  figure  are 
curves  calculated  using  an  extended  Donaldson  analysis 
reported  in  reference  (u) . For  those  values  of  Mach  number 
for  which  several  data  points  are  available  (M  « 5.0,  5.8, 
and  6.8),  the  results  indicate  that  increasing  values  of 
heat  transfei  parameter  have  little  effect  on  the  skin 
friction  ratio.  On  the  ether  hand.  Van  Driest  and  Monaghan 
(references  (s),  (y))  predict  an  increase  in  the  skin  friction 
ratio  up  to  about  10  percent  for  the  highest  heat  transfer 
case  Investigated  in  the  present  experiments.  It  should  be 
noted,  however,  that  it  is  implied  in  those  analyses  that 
the  edge  of  the  laminar  sublayer  occurs  at  fixed  values  of 
u"*"  and  y* , and  the  slope  of  the  turbulent  portion  of  the 
u"**,  y^  curve  remains  constant  regardless  of  Mach  numlor  or 
hoat  transfer  conditions.  As  pointed  out  previously,  the 
results  given  in  Table  I show  that  these  assumptions  are  not 
generally  valid.  This  is  more  evident  from  Figure  16  in 
which  the  values  of  Uj^  ■=  arc  plotted  as  a function  of 

the  heal  transfer  p.-iramc  tcr . Also  shown  are  the  associated 
theoretical  curves  derived  from  reference  (u).  These  figures 
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indicate  that  Mach  number  has  little  influence  on  the 
variation  of  Uj^  * with  (T^  - It  appears, 

however,  from  these  results  that  caution  should  be  exercised 
in  using  theories  which  embody  constant  values  of  u^  « y^ 

for  cases  of  heat  transfer. 

2S.  The  variation  of  with  Mach  number  is  shown  in 

Figure  17  together  with  the  diiect  skin  friction  measurements 
of  Coles  (reference  (z)}  and  the  deduced  skin  friction  data 
of  references  (a),  (q) , (r),  (t),  (y),  (z) , and  (aa).  In 
the  preparation  of  this  figure  the  data  for  M » 5 and  6.8 
were  reduced  to  the  case  of  zero  heat  transfer  by  fitting 
the  results  of  Figure  15  with  a straight  line  using  the 
method  of  least  squares.  For  the  other  Mach  numbers  the 
zero  heat  transfer  points  were  obtained  by  estimating  the 
slope  of  a straight  line  'hrough  the  data  consistent  with 
the  M c 5 and  6.8  results.  The  position  of  the  data 
points  on  this  figure  was  fixed  by  plotting  the  variation  of 
Cj/Cji  against  Re^  for  each  Mach  number  and  selecting 

the  values  of  at  a constant  value  of  Re^  of  8,000. 

This  value  of  Ro^  was  selected  because  it  represented  a 

Reynolds  number  for  which  the  greatest  overlap  of  data 
existed.  It  is  apparent  that  the  Cj/c^^  variation  with 

.Mach  number  is  a smooth  continuous  curve  for  a constant 
Reynolds  number.  Also  shown  in  this  figure  are  the  theoreti- 
cal predictions  of  relcrcnces  (q)  and  (u)  for  a constant 
Reynolds  number  of  8,000  and  zero  heat  transfer.  The  close 
fit  of  the  experimental  data  with  the  theoretical  curves 
indicates  that  the  theories  of  references  (q)  and  (u)  may  be 
used  to  predict  with  good  accuracy  the  variation  of  skin 
friction  coefficients  with  Mach  number  for  the  zero  heat 
transfer  case. 


CONCLUDING  REMARKS 

30.  Detailed  investigations  of  turbulent  boundary  layer 
velocity  and  temperature  profiles  have  been  made  on  a nozzle 
wall  of  tne  NOL  12  x 12  cm  Hypersonic  Tunnel  No.  4 at  Mach 
numbers  of  5.0,  5.8,  6.8,  7.7,  and  8.2  for  varying  rates  of 
surface  heat  transfer. 

31.  The  experimental  results  show  that  turbulent  boundary 
layer  profiles  in  hypersonic  flow  qualitatively  resemble 
turbulent  boundary  layer  profiles  in  incompressible  flow  in 
many  details. 

Si.  All  velocity  profiles  measured  can  be  fitted  with  a 
powrr  law  in  the  outer  turbulent  portion. 
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33.  As  in  the  case  of  incompressible  flow,  the  turbulent 
portion  of  the  profiles  differs  in  shape  from  the  logarithmic- 
law  velocity  profile.  The  discrepancy  increases  with 
increasing  heat  transfer. 

34.  The  dispersing  effects  of  Reynolds  nximber  and  heat 
transfer  parameter  exhibited  in  the  u*,  y*  velocity  profile 
representation  for  the  outer  turbulent  portion  could  bo 
removed  by  Lasing  the  coordinates  on  local  properties  at  each 
point  in  the  boundary  layer  and  dividing  the  distance  parame- 
ter by  the  respective  shape  parameter  (11) . The  data  for  any 
one  llach  number  then  fall  on  a single  curve.  The  velocity 
profile  curves  show  a systematic  displacement  with  Mach  number. 

35.  The  velocity  profiles  in  the  laminar  sublayer  are  linear. 

The  ratio  of  laminar  sublayer  thickness  to  total  boundary 
layer  thickness  decreases  slightly  with  increas- 

ing heat  transfer  rate  but  increases  considerably  with 
increasing  Mach  number  at  the  same  heat  transfer  rates. 

3G.  The  u*  = y”^  value  at  the  edge  of  the  laminar  sublayer 
for  the  zero  heat  transfer  case  is  close  to  the  incompressible 
flow  value.  However,  this  u'*'  « y'*’  value  is  found  to  increase 
with  increasing  heat  transfer  for  any  given  Mach  number. 

37.  The  ratio  of  the  thoimal  sublayer  thickness  to  the 

total  boundary  layer  thickness  was  found  larger 

than  iJ thermal  sublayer  tijickness  ratio  is 
unaffected  by  heat  tr.nnsfcr  but  somewhat  affected  by  Reynolds 
number . 

38.  The  data  strongly  support  the  applicability  of  Reynolds 
analogy  for  turbulent  I f^undary  layers  in  hypersonic  flow. 

30.  Values  of  skin  friction  coefficients  calculated  from 
the  heat  transfer  measurements  by  using  Reynolds  analogy 
agree  to  within  5 percent  with  those  skin  friction  values 
deduced  from  the  velocity  profile  slope  in  the  laminar 
sublayer . 

40.  For  the  test  range  of  the  heat  transfer  parameter 

(Tp  - T^)/Tp  (approximately  0 to  0.5),  values  of  the  skin 

friction  coefficient  are  found  to  be  only  slightly  affoct'Sd 
by  heat  transfer  rate  and  arc  in  accord  with  the  direct  .skin 
friction  measurements  of  Coles  and  other  investigators  when 
the  results  are  based  on  a single  value  of  Reynolds  number. 

41.  The  present  results  extend  the  range  of  available  skin 
friction  and  heat  transfer  data  to  a Mach  number  of  8.2. 

42.  The  detailed  measurements  of  velocity  and  temperature 
distributions  across  turbulent  bound.ary  layers  in  hypersonic 
flow  pi'csonted  here  enlarge  the  fund  of  available  data  for 
turbulent  boundary  layers  in  compressible  flows. 
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FIG.  i MACH  NUMBER  VARIATiG-N  ACROSS  BOUNDARY  LAYER  FOR  FREE- 
STREAM  MACH  NUMBERS  OF  5 AND  6.8  AND  VARIOUS  RAl  tS  OF 
HEAT  TRANSFER 


FIG.  6 VELOCITY  VARIATION  ACROSS  BOUNDARY  LAYER 
FOR  FREE-STEAM  MACH  NUMBERS  OF  5 AND 
6.8  AND  VARIOUS  RATES  OF  HEAT  TRANSFER 

VAl.l'ES  between  0<r<OS  BJtSEDON  INT£Rf>0».ATtn  E,  VALUES 

29 
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FIG.  9 SEMI-LOGARITHMIC  REPRESEN  I'AT ION  OF  BOUNDARY  LAYER 
VELOCITY  PROFILES  FOR  A FREE-STREAM  MACH  NUMBER 
OF  5.8  AND  THREE  RATES  OF  HEAT  TRANSFER 
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FIG  10  SEMI- LOGARITHMIC  REPRESENTATION  OF  BOUNDARY  LAYER 
VELOCITY  PROFILES  FOR  A FREE- STREAM  MACH  NUMBER 
OF  6 fi  AND  FOUR  RATES  OF  HEAT  TRANSFER 
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FIG.  12  SEMI-LOGARITHMIC  REPRESENTATION  OF  TURBULENT  PORTION 
OF  BOUNDARY  LAYER  PROFILES  IN  o*yVH  COORDINATES  FOR 
VARIOUS  MACH  NUMBERS 
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FIG.  16  VARIATION  OF  u\  = y\  VALUES 
WITH  HEAT  TRANSFER  PARAMETER 
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FIG.  17  VARIATION  OF  SKIN  FRICTION  COFFFICIENT  WITH  MACH 
NUMBER  FOR  ZERO  HEAT  TRANSFER  AND  A Rep  VALUE 
OF  8,000 

• THE  DATA  PLOTTED  ARE  I^ANSKIN  FRK,  I ION  KATIOS  ^F/Cpi 
»•  THE  DATA  PLOTTED  WERE  03TAif«c0  FROM  A STRA’CHT-LINE  EXTRAPOLATION. 
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